Iranian Journal of Numerical Analysis and Optimization 
Vol. 13, No. 4, 2023, pp 728-746 _— 


https: //doi.org/10.22067 /ijnao.2023.81484.1229 _ 

https://ijnao.um.ac.ir/ 

How to cite this article ro} 
BY NC Research Article 


Numerical study of sine-Gordon equations 
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Abstract 


The nonlinear space time dynamics have been discussed in terms of a hyper- 
bolic equation known as a sine-Gordon equation. The proposed equation 
has been discretized using the Bessel collocation method with Bessel poly- 
nomials as base functions. The proposed hyperbolic equation has been 
transformed into a system of parabolic equations using a continuously dif- 
ferentiable function. The system of equations involves one linear and the 
other nonlinear diffusion equation. The convergence of the present tech- 
nique has been discussed through absolute error, L2-norm, and L,o-norm. 
The numerical values obtained from the Bessel collocation method have 
been compared with the values already given in the literature. The present 
technique has been applied to different problems to check its applicability. 
Numerical values obtained from the Bessel collocation method have been 
presented in tabular as well as in graphical form. 
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1 Introduction 


Nonlinear partial differential equations have wide applications in different 
branches of science and engineering, which helps to understand the diversity 
of physical phenomena in a logical manner. The nonlinear wave equations 
such as Klein—Gordon and sine-Gordon equations have wide applications in 
optics, plasma physics, quantum mechanics and fluid mechanics, and so on. 
Finding the analytic solution to these problems is a tedious task due to 
the complexity of the nonlinear terms. The numerical approximation of the 
solution in terms of the discrete set of points is often desirable by researchers 
due to the simplicity of the numerical computations. 

The derivation of the Klein-Gordon equation is a generalization of the 
Schrédinger equation. It was named after the physicists Oskar Klein and Wal- 
ter Gordon. They together in 1926, proposed that relativistic electrons can 
be described by the Klein—Gordon equation during the research for the equa- 
tion describing de Broglie waves. Schrédinger considered the Klein—Gordon 
equation as a quantum wave equation [18, 19]. Klein—Gordon equation plays 
a significant role in many scientific applications, such as nonlinear optics and 
quantum field theory, and solid state physics. Klein—Gordon equation in one 
dimension can be considered as 

2 2 
Te Page + hu) = fl6.0) (1 


where (£,t) € (4,6) x (0,7), fi(y) is nonlinear force, and / is a constant. 
The sine-Gordon equation is a special case of (1) for fi(y) = sin(y) and 
fo(€,t) = 0. Interestingly, the sine-Gordon equation was discovered sepa- 
rately in 1939 by Frenkel and Kontorova while studying the propagation of 
slip in an infinite chain of elastically bound atoms lying over a fixed chain. 
Later it was found that sine-Gordon is a special case of Klein—Gordon equa- 
tion. The sine-Gordon equation in one dimension can be described as 


ae ae em) =O (2) 
and the initial and boundary conditions can be described as 


y(E, 0) = o1(€) and y(E, 0) = a(€), 
y(Ea,t) = G3(t) and y(f,t) = da(t). 


The present problem of the equation is a continuum model for waves in me- 
chanical systems, coupled-pendulum, the study of the domain wall dynamics 
in magnetic crystals, magnetic-flux propagation in large Josephson junctions, 
propagation of crystal dislocations in solids, propagation of ultra-short op- 
tical pulses in optical fibers, as a nonlinear effective field theory for strong 
interactions in particle physics, and so on; see [8, 16, 18, 19, 28, 37]. 


(3) 
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A variety of numerical techniques have been developed by different re- 
searchers to study the behavior of nonlinear sine-Gordon equation, such 
as finite difference method, inverse scattering method, auxiliary equation 
method, spectral method, pseudo-spectral method, tanh-sech method, Ado- 
mian decomposition method, sine-cosine method, Jacobi elliptic functions, 
Backlund transformation, Riccati equation expansion method, homotopy per- 
turbation method, and variational iteration method [15, 18, 33, 32, 40]. 

Method of characteristics and a leapfrog finite difference scheme [3, 24] 
were the first two methods developed to obtain the numerical solution of 
sine-Gordon equation. Strauss and Vazquez [36] developed a leapfrog finite 
difference, an implicit, energy-conserving scheme for the Klein—Gordon equa- 
tion. 

Apart from these, other numerical methods have also been developed for 
the solution of sine-Gordon equations, which include pseudospectral meth- 
ods and spectral methods. Pseudospectral methods include the split-step 
Fourier scheme [1, 2, 41] and spectral methods include energy-conserving, 
wavelet spectral method, Fourier scheme, Legendre spectral element method, 
and multiresolution analysis method based on Legendre wavelets [27]. Finite 
element methods is based on a collocation scheme using Legendre—Gauss— 
Lobatto points [29], cubic B-splines [31] and Petrov-Galerkin scheme [5]. 

In the present study, the Bessel collocation method (BCM) has been fol- 
lowed to study the behavior of the nonlinear sine-Gordon equation. The 
Bessel polynomials of degree n have been taken as base polynomials. To 
discretize the time direction, the temporal variable has been split by the 
introduction of a continuously differentiable function. It converts the wave 
equation into a system of equations involving one linear and the other non- 
linear equation. 

The present manuscript has been divided into six sections starting from 
introduction. The BCM has been described in section 2. The explanation of 
collocation points as well as the implementation of the collocation technique 
has been described in sections 3 and 4, respectively. Convergence analysis 
has been discussed in section 5, whereas the numerical application has been 
given in section 6. 


2 Bessel collocation method (BCM) 


The collocation method belongs to the general class of approximate methods, 
known as weighted residual methods. In this method, the residual is set 
orthogonal to the weight function. In an orthogonal collocation, the trial 
function y(€,t) is represented in a series of known polynomials with unknown 
coefficients [13, 22, 38]. The residual is set equal to zero at the collocation 
points. 

On the basis of the implementation of the trial function, the collocation 
technique can be classified into three categories. If the trial function satisfies 
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the differential equation €V(y) = 0 with volume V, then it is termed as inte- 
rior collocation. If the trial function satisfies the boundary @B(y) = 0, where 
B is the boundary adjoining volume V, then it is termed as boundary collo- 
cation. If the trail function satisfies neither the equation nor the boundary 
and is adjusted to both, then it is termed as a mixed collocation. 

The choice of base function is the first important step in the technique of 

collocation. In the present study, Bessel polynomials of order n have been 
chosen as a trial function, and the technique is called the BCM. 
During the study of problems in dynamic astronomy to solve the Kepler’s 
problem, a German astronomer Bessel in 1824, introduced Bessel polyno- 
mials, which are the solution to a second order boundary value problem. 
These polynomials can be written in terms of limit confluent hypergeomet- 
ric function 9 F,. The details of these hypergeometric functions are given in 
(7, 25, 34, 39]: 


n 


Inlé) = se goFi(-in + 1; 42). (4) 
The Bessel coefficients also follow from the power series expansion for small 
values of € 
pare 1; —72) =, 
lim €-"Jn(€) = =, 
£30 2"n! 


which shows that as € — 0, the Bessel coefficient J,,(€) approaches to xT 
The first order derivative of the Bessel function is defined as 


d n en 

Pra In (E)) = €" In_1(8), 
Dees, a 
ae Jn(€)) = —E-" Ing (€) 


3 Collocation points 


The next step is the choice of collocation points. It is an important part of 
the collocation technique. In this study, instead of taking the uniform points, 
the zeros of orthogonal polynomials, such as Jacobi polynomials, have been 
taken as collocation points. Legendre and Chebyshev polynomials are special 
cases of Jacobi polynomials, and the zeros of these orthogonal polynomials are 
preferably taken as collocation points. Runge’s divergence formula also states 
that nonuniform collocation points give less error as compared to uniform 
collocation points. 


Theorem 1. [26] If Q,.(€) form a simple set of real polynomials and w(€) > 0 
ona< & <b, then the necessary and sufficient condition that the set Q,(§) 
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is orthogonal with respect to the w(€) over the interval a < € < b is that 


b 
: w(f)a* QO, (€)dé = 0, k=0,1,2,3,...,(n—1). 


Theorem 2. [26] If the simple set of real polynomials Q,,(€) is orthogonal 
with respect to the weight function w(€) > 0 on the interval a < € < b, then 
the zeros of Q,,(€) are distinct and lie in the interval a < € < b. 


Since Q,,(€) is a polynomial of degree n, then it has exactly n roots, 
multiplicity counted, such that the roots are distinct and all lieina<€ <b. 

Usually, the collocation points are selected from the Legendre or Cheby- 
shev polynomials, and these polynomials are also particular cases of Ja- 
cobi polynomials. The details of the collocation points are given elsewhere 
[4, 6, 11, 12, 14, 17, 20, 35]. The zeros of Chebyshev polynomials have been 
taken as collocation points: 

g = cos (HUY), 7 =1,2,...,n41. 

The Chebyshev collocation points have been transformed from the interval 
[—1, 1] to [&.,] using one to one correspondence. 


4 Implementation of BCM 


To discretize the given problem, BCM is applied in space direction. To 
apply BCM, (2) has been split in time direction. For this purpose, a new 
continuously differentiable function z(€, t), differentiable with respect to t has 
been introduced: 


oe = 26,1), (E,t) € (€a,€) x (0,7), 


ee (5) 
7: = - —sin(y(€,t)),  (€ t) © (€a,6) x (0,T). 


To apply BCM on a system of equations defined by (5), the two functions 
y(€,t) and z(€,t) have been approximated in terms of Bessel polynomials as 


n+1 
y(E,t) = », Ji(E)ex(t), 

a (6) 
2(€,t) = >i Ji(E)di(t), 


where J;(€) are ith order Bessel polynomials. To simplify (6), the Bessel 
polynomials can be rewritten as suggested by [10, 42, 43, 44, 45, 46, 47, 48]: 
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n+1 


wlét) = re Rei(t), 


n+1 


z(é,t) = » é'' Rd; (t), 


(7) 


where R is a square matrix of order (n+ 1) x (n+ 1) and c;(t) and d;(t) are 


the unknown coefficients of t, which are to be determined. 
For n being an odd integer, R is defined as 


1 
010120 ; 0 0) 0 
o oltal 2 o 2 
= 1. 
11122 ol2!22 7 
O 112123 . 013123 e 
212104 0 113124 z o 
R= 
n—-1 n—-3 
(-1) 2 6 (-1) 2 0 Sane ee 
(Beer (222 en O!(n—1)!2%— 
n-1 n—3 
(-1) 2 (-1) 2 
0 > 0 ae 0 
(AE (42 EB an 
However, for n being an even integer, the matrix R can be written as 
1 
dio120 (0) ie) 0 0 0 
0 onrat ) 0 0) 0 
— 0 a 0 0 oO 
1!1!2 1 o!2!2 < 
e 112123 013123 9 o 
solor 0 aii 0 0) 0 
21212 1!3!2 
R= 
‘n-2 ‘ ‘nA 
(-1) 2 F (»2 1 
HEIR PRAT = + = aera = 
Apert (23 P int Oman 
n Mee 
(-1)2 ‘ (-1) 2 F F 1 
L ye ye (2278 en a OmIa7 | 


At jth collocation point, (7) can be written as 


n+1 
yi, => Ghat), j=1,2,...,0+ 1, 
i=1 


n+1 
AG ith= > & Rat), 7=1,2,..,0+1. 
i=1 


Now, rewrite (8) in a matrix form at jth collocation point 


w=), [zs] = XM[d I, (9) 
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where X = [7 |R and y; represents the value of y at jth collocation point. 
We have 
[X]"[ys] =[el,  [X)* [ze] = Id. (10) 


Substituting collocation coefficients from (10) in (7) results in 


n+1 ; 

y(€, t) = S- ‘ame ue 
a 

2(€,t) = ee ae 


i=l 


(11) 


The first and second order derivatives of y(&,t) with respect to € can be 
obtained as 


Oy = . i-2y-1 
= (12) 
ay n+1 


ae = Soli — 1) — 2)e°-F XT. 


i=l 


Using the discretized forms of y(€,t) and z(€,t) in (5) leads to the following 
system of equations: 


de n+1 (13) 


A = BY~ Buy: — sin(y;), 
1=1. 


where j = 2,3,...,n. 

In the above coupled form of equations, Aj; and B,,; are first and second order 
discretized forms of derivatives of y(€) with respect to € at jth collocation 
point, respectively. Boundary conditions for both y(€,t) and z(&, ¢) configura- 
tions assumed to be y(a,t) = y1 = Ya, y(b,t) = Yn41 = Yo, 2(a,t) = 21 = Za, 
and z(b,t) = Zn41 = 2. The matrix representation of (13) of sine-Gordon 


can be written as: ie 
See ITO\lZ 
[2] =[o8][¥]-F as 


dt 
In the above system of equations, there are n — 1 collocation equations and 
two boundary conditions for each function y(€,t) and z(€,t), respectively. It 
results in 2(n—1) collocation equations in total and four boundary conditions. 
There is no effect of boundary conditions in the matrix representation of (14) 
as they are in scalar form and get merged into F’. Moreover, O represents 
the zero matrix, and I represents the identity matrix in relation (14). Also, 
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000...0 111...1 
000...0 111...1 

=l]....°. ; I=]...... 

(n—1)x(n—-1) 


000...0 11i1...1 


(n—1)x(n—1) 
Y2 22 
¥3 23 
Y= , Z=|. ; 
Yn (n-1)x1 en (n-1)x1 
B(2,2) B(2,3) B(2,4) ... B(2,n) 
. B(3,n) 


and 


sin(Yn) (n-1)x1 


The matrix representation corresponds to a nonlinear system of 2(n — 1) 
equations form block matrix structure (14). The left-hand side includes the 
column vector of time derivatives of functions y and z, respectively, and the 
right-hand side does not include any product with the inverse of a coefficient 
matrix, which helps to reduce the stiffness of the system of equations. The 
system of ordinary differential equations is solved using MATLAB with the 
ode23s module. 


5 Convergence analysis 


Theorem 3. [26] If {Q,(€)} represents for simple set of polynomials and if 
Y(E) is a polynomial of degree m, then there exist constants a, such that 


VE) = So ax Qe (6). 
k=0 
The a,’s are functions of & and any parameter involved in Y(E). 
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Theorem 4. [25]There exists a unique polynomial P,,(€) of degree n, which 
assumes prescribed values at n+ 1 distinct points £) < £; <--- < &. 


Theorem 5. [25] Given any interval a < € < 6, real number ¢ > 0 and any 
real valued continuous function f(€) on a < € < b, there exists a polynomial 
P(€) such that 


Il F(€) — P(€) I<. 


To study the convergence behavior of orthogonal collocation scheme, the 
definition given by [21] is quoted here: 
Consider a family of mathematical problems parametrized by singular pertur- 
bation parameter e, where « lies in the semiopen interval 0 < ¢ <1. Assume 
that each problem in the family has a unique solution denoted by y-, and that 
each y~- is approximated by a sequence of numerical solutions {(Ye, ou) eee 
where Y- is defined on the QN representing the set of points in R, and N 
is the discretization parameter. Then the numerical solutions Y- are said to 
converge to the exact solution y-, if there exists a positive integer No and 
positive numbers G and p, where No, G, and p are all independent of N and 
€, such that for all N > No 


supo<e<i|Ye — ie | <GN-?, 


Here p is the rate of convergence and G is the error constant. It shows that 
the rate of convergence in the case of the collocation technique depends upon 
the number of collocation points. 

Let y(&,t) be the exact solution, and let y,(&,t) be the approximate so- 
lution. The absolute error between the exact and approximate solution is 
calculated as 

Ea =| y(&,t) — ynl§,t) |. (15) 


The error in terms of Lo-norm and L,.-norm has been calculated with re- 
spect to the weight function w(&) such that 


n+1 


lly — yalls = S— lwilE)(y — yndi I, (16) 


i=1 
where y(&,t) represent analytic solutions and y;,(€,t) represent approximate 
solutions [22]. The error between exact and numerical values has been shown 


by e=Y—-—Yn- 
L2-norm is said to converge to the exact solution if ||y — y,||2 > 0 as 
n — oo. Thus 


llell2 =Ily — valle, (17) 


Similarly, in £..-norm for ||y — yp||, it has been taken as 
lly — Ynlloo = max |(y— yn)al, 1=1,2,3,...,.n+1, (18) 


llelloo = Il¥ — Yalloo- (19) 
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6 Numerical examples 


To verify the applicability of BCM, the scheme has been applied on different 
nonlinear hyperbolic equations. 


Example 1. Consider the sine-Gordon nonlinear hyperbolic equation from 
the generalized (2) coupled form of two interacting configurations y(&,t) and 


Sy (é,t) = 2(€,t): 


Ory Ory 
Ot 0€2 sin(y), 
by defining 
oy 
Ot , 
Oz Ory 


with respect to the initial conditions 


y(,0) = 4tan~*(exp(g€)), 


exp(g&) 
2(€,0) = —4eg (1 + exp(2g€)) ’ 


Then exact solutions of the above equations have been obtained as 


y(€,t) = 4tan™*(exp(g(€ — ct))), 


where € € [—3,3], g = aa and c = 0.5, and boundary conditions can be 


extracted from the exact solution [9, 23, 30]. 

The values of y(€,t) are calculated for different numbers of collocation 
points and compared with the exact solution in Table 1. The numerical 
values have been presented for a fixed value of € = 0 and different values 
of time in Table 1. The numerical results are found to be enough close to 
the exact solution. It is also observed that no particular change occurs in 
numerical values after 25 collocation points. A graphical representation of 
experimental results with respect to the time and collocation points has been 
presented in Figure 1. 

A comparison of error in terms of L2-norm and L..-norm at different 
numbers of collocation points has also been given in Table 2 for different 
values of time. 

A graphical representation of error in terms of Z2-norm and L.-norm 
at 25 and 31 collocation points with respect to the time values have been 
presented in Figures 2 and 3, respectively. It can be analyzed from Figures 
2 and 3 that error decreases with the increase in collocation points. From 
these figures, it can also be analyzed that the decrease in error is not so large 
and almost similar at 25 and 31 collocation points. 
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Numerical results have also been compared with [9, 23, 30] and have been 
discussed in Table 3. It is observed that the results obtained by BCM are 
better and give less error. 


Example 2. The sine-Gordon nonlinear hyperbolic equation from the gener- 
alized (2) in coupled form of two interacting configurations y(€,t) and z(€,t) 
with respect to the initial conditions 


y(é,0) = 0, 
2(€,0) = 4sech(&), 
the exact solution of equation has been obtained as 
y(€,t) = 4tan—!(t.sech(€)), 


sech(€) 
1+ t?sech?(€)° 


The boundary conditions have been taken from the exact solutions [9, 23, 30]. 


2(€,t) =4 


The numerical values calculated at 17 and 19 collocation points have been 
compared with the exact values and are presented in Table 4 for fixed € = 0.5 
but at different values of the time. It has been observed from Table 4 that the 
numerical results are close enough to the exact solutions. It is also observed 
that no particular change in numerical values occurs after 17 collocation 
points. The graphical representation of numerical values with respect to the 
time and collocation points has been presented in Figure 4. 

A comparison of error in terms of Lo-norm and L,.-norms at 17 and 
19 collocation points has been given in Table 5 at different time levels. The 
graphical representation of error in terms of D2-norm and L..-norm at 17 and 
19 collocation points with respect to the time has been presented in Figures 
5 and 6, respectively. It is observed from Figures 5 and 6 that at 17 and 19 
collocation points, it does not make much difference in numerical results, but 
if still count the difference at 19 collocation points, it shows better results 
than 17 collocation points. 

A comparison of numerical results with [9, 23, 30] has been discussed in 
Table 6 and found to be close enough to be accepted. 


7 Conclusion 


The given nonlinear sine-Gordon equation has been solved successfully by us- 
ing the BCM over Chebeshev collocation points. By the above analysis, the 
proposed method of BCM is proved to have some desired and popular fea- 
tures, such as high order accuracy and preserving energy conservation. Con- 
sistency and convergence of the computational technique have been obtained 
by computing the results of numerical solutions with analytic solutions. Error 
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analysis in terms of L2- and L..-norms with respect to the weight function 
employed showed that the Bessel collocation approach is very stable, and the 
results obtained by this approach are consistent and convergent. 


Table 1: Comparison of absolute error (Eq) of Example 1 at € = 0 for different numbers 
of collocation points 


t Ea at8 Eqgat1l6 Egat25 Ea at 31 
collocation collocation collocation collocation 
points points points points 

0.1  1.2362e-02 3.6988e-04 1.3449e-07 4.0068e-07 
0.2  2.4923e-02 7.4329e-04 9.8035e-07 1.9933e-07 
0.3 3.7862e-02 1.1126e-03 2.8106e-06 1.9219e-06 
0.4 5.1315e-02 1.4626e-03 5.2349e-06 1.9219e-06 
0.5 6.5366e-02 1.7749e-03 7.3092e-06 5.0514e-06 
0.6 8.0030e-02 2.0225e-03 7.9431e-06 1.6903e-06 
0.7  9.5257e-02 2.1711e-03 6.4173e-06 1.3360e-06 
0.8 1.1093e-01 2.1925e-03 2.7694e-06 5.4602e-06 
0.9  1.2689e-01 2.0728e-03 2.1323e-06 5.0533e-06 
1.0  1.4291e-01 1.8073e-03 6.8670e-06 6.3394e-06 


Table 2: Comparison of error for y(€,t) of Example 1 at different collocation points 


t At 25 At 31 
collocation collocation 
points points 
lelloo llell2 llelloo llell2 


0.1]6.7963e-07 4.3658e-07 |1.2869e-06 0.0000e-00 
0.2}2.4566e-06 1.2844e-06 |3.5131e-06  1.2717e-07 
0.3}4.7272e-06  1.9502e-06 |4.9638e-06 1.1174e-07 
0.4/6.7848e-06  2.5348e-06 |6.4823e-06 1.6430e-07 
0.5)7.3091e-06 3.7359e-06 |7.1190e-06 2.9715e-07 
0.6] 7.9430e-06 5.3058e-06 |7.7593e-06 2.8624e-07 
0.7|6.8875e-06 6.2351e-06 |7.9344e-06  1.3433e-07 
0.8] 7.9192e-06 5.8165e-06 |8.1705e-06 0.0000e-00 
0.9}9.8062e-06 4.2143e-06 |8.1533e-06 2.8093e-07 


Table 3: Comparison of ||e||2 and |{e||.. calculated by Bessel collocation for y(€,t) with 
different techniques of Example 1 


t |Dehgan & Shokri [9]| Mittal & Bhatia [23]|Shukla & Tamsir [30]| Bessel collocation 
llellos llelle Ilelloo llell2 llelloo llell2 llelloo__ lela 

0.25 /4.95e-06 1.76e-05 |4.90e-05 3.66¢e-05 |9.61le-06 5.67e-06 |4.44e-06 1.46e-07 
0.50 /8.42e-06 4.31le-05 |7.55e-05 9.00e-05 |1.10e-05 8.39e-06 |7.38e-06 2.97e-07 
0.75 |1.65e-05 8.25e-05 |1.43e-04 1.60e-04 |1.26e-05 1.05e-05 |7.99e-06 0.00e-00 
1.00 |2.51e-05 1.27e-04 |2.10e-04 2.27¢e-04 |1.44e-05 1.24e-05 |1.49e-05 2.81e-07 
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Table 4: Comparison of absolute error (Eq) of Example 2 at € = 0.5 for different numbers 
of collocation points. 


t Ea 


at 8 


Ea,at16 E,at17 EF, at 19 


collocation collocation collocation collocation 


Pp 


oints points 


points 


points 


0.1 4.6603e-06 1.0754e-10 4.4239e-10 5.3072e-11 
0.2 2.9466e-05 8.3158e-10 6.1702e-11 2.3970e-10 
0.3 6.5426e-05 2.2862e-09 1.1462e-09 6.9454e-10 
0.4 7.9905e-05 3.3772e-09 2.5359e-09 1.1188e-09 
0.5 5.2043e-05 3.7229e-09 3.3682e-09 1.4447e-09 
0.6 4.2174e-05 3.7595e-09 3.1963e-09 1.4496e-09 
0.7 3.7627e-05 3.0807e-09 2.0244e-09 1.0735e-09 
0.8 3.1964e-05 7.4804e-10 3.7364e-10 4.0962e-10 
0.9 2.4734e-05 3.0871e-09 1.3676e-09 4.7658e-10 
1.0 3.1170e-05 6.9712e-09 3.3311e-09 1.5792e-09 


Table 5: Comparison of error for y(€,t) of Example 2 at different collocation points 


At 17 
collocation 
points 


At 19 
collocation 
points 


lle| 


loo llell2 


llelloo llell2 


3.6956e-10 
7.5355e-10 
2.5508e-09 
4.2276e-09 
4.5143e-09 
3.0879e-09 
2.2236e-09 
2.2897e-09 
4.4816e-09 


2.0235e-10 
4.0045e-10 
1.6149e-09 
2.6935e-09 
3.0129e-09 
2.4971e-09 
1.4627e-09 
1.1368¢-09 
2.4093e-09 


1.2951e-10 
4.6931e-10 
1.3351e-09 
1.9375e-09 
2.0691e-09 
1.5904e-09 
1.2266¢e-09 
8.0804e-10 
1.5020e-09 


7.3902e-11 
2.9364e-10 
8.3540e-10 
1.2585e-09 
1.4681e-09 
1.3228¢-09 
9.1829e-10 
5.1226e-10 
7.6839e-10 


Table 6: Comparison of ||e 
different techniques of Exa: 


J2 and |le||oo calcula: 
mple 2 


ed by Bessel collocation for y(€,¢) with 


Dehgan& Shokri [9] 


Mittal& Bhatia [23] 


Shukla& Tamsir [30] 


Bessel collocation 


llelloo llell llelloo llell Ilelloo llell2 llelloo _lella 
0.25 ]5.89e-06 3.91e-05 [2.32e-05 1.18e-05 5.46e-06  2.43e-06 ]9.06e-10 5.64e-10 
0.50 [2.01e-05 1.30e-04 [4.11e-05 4.19e-05 |7.39e-06 5.54e-06 |1.88e-09  1.44e-09 
0.75 |3.63e-05 2.35e-04 |1.02e-04 7.78e-05 |7.78e-06 6.45e-06 |1.02e-09 6.87e-10 
1.00 [5.07e-05 3.27e-04 [1.64e-04 1.30e-04 | 8.75e-06 _7.84e-06 |2.54e-09 _1.53e-09 
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Figure 1: Graphical representation of y(£,t) of Example 1 
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|| — © —llell, at 31 collocation points 


Figure 2: Graphical representation of comparison of error in form of L2-norm with 
respect to time and number of collocation points of Example 1 
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Figure 3: Graphical representation of comparison of error in form of Lo.-norm with 
respect to time and number of collocation points of Example 1 
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Figure 4: Graphical representation of y(€,t) of Example 2 


S107? 


——*— llellg at 17 collocation points 
4) |— © —llellz at 19 collocation points 


fo) 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 
Time 


Figure 5: Graphical representation of comparison of error in form of L2-norm with 
respect to time and number of collocation points of Example 2 
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Figure 6: Graphical representation of comparison of error in form of Loo-norm with 
respect to time and number of collocation points of Example 2 
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